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Abstract 


Many of the new physics models predicts a light Higgs boson similar to the Higgs boson of 
the Standard Model (SM) and also extra scalar bosons. Beyond the search channels for a SM 
Higgs boson, the future collider experiments will explore additional channels that are specific to 
extended Higgs sectors. We study the charged Higgs boson production within the framework of 
two Higgs doublet models (THDM) in the proton-proton collisions at the FCC-hh collider. With 
an integrated luminosity of 500 fb _1 at very high energy frontier, we obtain a significant coverage of 
the parameter space and distinguish the charged Higgs-top-bottom interaction within the THDM 
or other new physics models with charged Higgs boson mass up to 1 TeV. 
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I. INTRODUCTION 


The Higgs sector of the Standard Model (SM) is in the minimal form, which contains 
one complex isospin doublet of scalar fields (0 + ,0°), resulting in one physical neutral C 
even Higgs boson (h°) which has been discovered at the LHC by ATLAS [1] and CMS 
collaborations. However, there are many possibilities for the extension of the Higgs sector, 
introducing further multiplets of scalar fields, which might be singlets, doublets, or triplets 
of the symmetry groups. The extended Higgs sectors are more relevant to the neutrino 
mass, baryogenesis and dark matter. It is thus natural to consider scenarios with additional 
complex scalars such as two Higgs doublet model (THDM) jjij]. The scalar fields ((t>\ and 


02 ) couplings to the up-type, down-type and c 
identified for discrete types of the THDM 4, 5 


rarged lepton SU(2)l singlet fermions can be 
. Model of type-I are the one in which all SM 


fermions couple to a single scalar field ( 02 )- In type-II model down-type quarks and charged 
leptons couple to one scalar field (0i), while the up-type quarks and neutrinos couple to 
the other scalar field (02). Furthermore, in the model of type-III the quarks couple to one 
of the scalar field (0i), while leptons couple to the other (0 2 ). In the model of type-IV, 
the couplings of scalar (0 2 ) to up-type quarks and charged leptons, and the couplings of 
scalar (0x) to down-type quarks and neutrinos are present. Here, we consider two types of 
THDM-I and THDM-II for type-I and type-II of two Higgs doublet model, respectively. 


The two Higgs doublets carry opposite hypercharges, the (0 1 and 0 2 ) scalar potential will 
contain the mixing parameters related to the mass. In this case, the Higgs doublets will have 
different vacuum expectations values (iq and rq). The massive SM gauge bosons acquire 
their masses from the expressions with the vacuum expectation value v = (vf An!) 1 / 2 . After 
the spontaneous symmetry breaking there appears five physical scalar particles: two neutral 
CP-even bosons h° and f/°, one neutral CP-odd boson A 0 , and two charged bosons H ± . 
Phenomenologically, the two Higgs doublet model includes the free parameters: mixing 
angle, the ratio of the vacuum expectation values (tan/3 = tq/tq), the masses of Higgs 
bosons. In the extended models of multiple neutral scalar bosons, the mixing between them 
would make it difficult to identify their properties. Therefore, it is important to study the 
charged scalar bosons, which could provide unique signatures to distinguish the models with 
extended Higgs sector. 


Constraints on the charged Higgs bosons in the THDM are given from both low energy 
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flavour experiments and high energy collider experiments. The direct lower bounds on the 
charged Higgs boson mass m#± come from LEP experiments. They were sensitive to the 
masses of charged Higgs boson up to about 90 GeV, in two decay channels of H + —» r + v 
and H + — > cs , and the exclusion limit on the mass independent of the admixture of these 
branching fractions was 78.9 GeV 4], Assuming the THDM-II model the lower bounds 
are 87.8 GeV for large tan/3 values using the decay mode H + —> t + v and 80.4 GeV for 
relatively low tan/3 values. Using the bounds m^o > 92 GeV and the characteristic relation 
m H ± = {m 2 A o + mw±y / 2 (as in the MSSM), one obtains the bound rn h± > 122 GeV. 
In the model THDM-II the mass of charged Higgs boson is constrained by the precision 
measurements of the radiative decay of B X s ^ by the low energy experiments. The 
lower bound on the mass m#± > 295 GeV are given for the model THDM-II. The decay 
B —>■ tv can also be used to constrain the charged Higgs parameters, being sensitive to 
(tan/3/m#±) 2 , which yields a lower bound m#± > 300 GeV for tan/3 > 40. A study of 
data analysis based on b —> sy [6] have excluded a mass range up to 380 GeV in a variety 
of interesting processes and BSM scenarios. However, these bounds can be relaxed if the 
MSSM or other new physics models contributes through the loop diagrams. 

Beyond the search channels for a Standard Model Higgs boson, the LHC experiments are 
exploring additional channels that are specific to extended Higgs bosons. ATLAS [7] and 
CMS js] collaborations have already performed a number of extended Higgs searches which 
exclude tan/3 > 50 in the range of heavy charged Higgs, m#± > 200 GeV. Recently, the 
CMS collaboration js] have put 95% C.L. exclusion limit on the mass of the charged Higgs 
boson in 180-600 GeV mass range. This search is performed at a center of mass energy of 
8 TeV with 19.7 fb _1 of data from pp — > t{b)H + and pp —> t{b)H~ production processes 
with H + — >■ r + v T decay mode. The ATLAS collaboration 10] have searched for charged 
Higgs bosons decaying through H ± —> process using the proton-proton collision data 

at y/s = 8 TeV with 19.5 fb 1 , the results exclude a large range of tan/3 values for charged 
Higgs boson masses in the range 80 —160 GeV, and exclude parameter space with high tan /3 
for the range of mass m#± = 200 — 250 GeV. 

The searches of the heavy Higgs bosons of the THDM have special challenge at present 
high energy colliders. One of the future international projects currently under consideration 
is the Future Circular Collider (FCC) [ll] which has the potential to search for a wide 
parameter range of new physics. The FCC-hli collider is to provide proton-proton collisions 
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at nearly an order of magnitude higher energy than the LHC, t 
energy is 100 TeV and the peak luminosity is 5 x 10 34 cm _ 2 s _1 

In this work, we study the processes pp —> ttb(b ) + A" for the signal and background 
at the FCC-hh collider. Our analysis is focused on the production of a pair of top quarks 
and associated bottom quark for the charged Higgs boson search within the THDM-I and 
THDM-II in the pp collisions at very high energy frontier. We have obtained the significant 
coverage of the parameter space at large integrated luminosity projections for the FCC-hh 
collider. We define the relevant expressions for the H~qq and H~l + u couplings as well as 
scalar-vector-scalar (H~W + h °, H~W + H°, H~W + A °) couplings in section II. We present 
the calculation for the decay widths and branchings ratios of the charged Higgs boson in 
section III. In section IV, we plot the production cross section according to the mass of 
charged Higgs boson for two types of THDM at the center of mass energy y/s = 100 TeV 
of the pp collisions. Finally, the results from the analysis of the signal and background are 
given in section V. 


le proposed centre-of-mass 

Q. 


II. CHARGED HIGGS BOSON COUPLINGS 

The magnitude of the couplings for H~fifj interactions are given by 

dH-qiqj = g\Vq iqj \[m qi cot 0(1 + 75 ) + T x m qj tan0(1 - 'Ts)]/(2v / 2 m w ) (1) 

9H-uii+ = glUuii^Txm^ tan 0(1 - 75 )]/( 2 V2m w ) ( 2 ) 

where V qiQj and U Vi i j are the CKM matrix elements in quark sector and PMNS matrix 
elements in lepton sector, respectively. The g is the weak coupling constant, the tan/3 
is the ratio of vacuum expectation vales of the Higgs doublets. The rn q and mi are the 
corresponding quark and lepton masses, respectively. We use the parameter T x to identify 
the type of THDM such that Tj = —cot0/tan0 denotes THDM-I and Tu = 1 denotes 
THDM-II. 

The couplings for H~W + h°, H~W + H° and H~W + A° interactions can be written by 

gH-w+ho = gcos(/3 - a)(p H - +p h o)/2 (3) 
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9h-w+ho = gsm(0 - a){p H - + p H o)/2 


( 4 ) 


Qh-w+ao = ig{PH- +Pao )/2 ( 5 ) 

where the p h o, p H o, and p H - are the four-momenta for neutral Higgs bosons and charged 

Higgs boson, respectively. We use the expression for the angle factor of cos 2 (f3 — a) = 
[(m^o/m|)(m^o/m|-l)]/[(m^o/m|-m^o/m|)(mio/m|+m^ 0 /m|-l)]. By implementing 


the relevant interaction vertices into the CalcHEP 


13] package, we calculate the decay width 


and branching ratios as well as the process cross section for different parametrizations (PI 
for rriA o = 100 GeV and m H o = m H - , PII for m^o = m#o = m#-) within the THDM-I and 
THDM-II. 


III. DECAY WIDTH AND BRANCHING RATIOS 

The partial decay widths of charged Higgs boson into fermionic channels can be calculated 
as 


P (H -A qiQj) = 


3<? 2 A 1 / 2 (m 2 H -, m 2 ., m 2 .) 


32ir7n‘y V ni 3 H - 


\V< 


<UQj I 


X 


(m 2 H - - m], - m 2 q .){T 2 m 2 q . tan 2 j3 + m 2 cot 2 P) - 4 T x m 2 q .m 2 q . 




( 6 ) 


<? 2 A 1//2 (mll_, m?., 0) 

T(H~ Vilj) = J f ’ 3 J ’ V ^hl 2 


327T m^rriff- 


Tl{m 2 H - - m\)m 2 h tan 2 ^ 


( 7 ) 


where A 1//2 is a kinematic factor of mass squared dimension and it can be defined as 
A(m 2 , m|, mf) = [(m 2 + m| — m§) 2 — 4m 2 m 2 ]. For the channels H~ —> W~h °, H~ —> W~H° 
and H~ —y W~A° we calculate the decay widths 

g 2 \ 1/2 {m 2 H -,mw,m 2 h0 ) cos 2 (P - a) 


T{H~ -A r h u ) = 


64:7im 3 H - 


x 


- 2{m 2 H - +m 2 0 ) + 


i m H- 


m 


2 N 2 


m 


IV 


( 8 ) 


t(h~ ->• rn°) = 


g 2 \ 1/2 (m 2 H -,mw,m 2 H o) sin 2 (/3 - a) 


647rm^_ 
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Figure 1. Decay width of H± boson depending on its mass for different values of tan/3 and 
parameter set PI of THDM-I. 
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Y{H~ -f W~A°) = 


fX 1 / 2 (m 2 H ,m 2 Wl m 2 A0 ) 
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( 10 ) 


The decay widths of the charged Higgs boson for the models THDM-I and THDM-II are 
presented in Fig. [TJHJ It is shown that the decay width for model THDM-I is larger than 
that for THDM-II in the considered parameter space. The decay width increases with the 
charged Higgs boson mass for parameter PI of THDM-I, while it slightly changes depending 
on the mass m#- > 300 GeV for parameter PII of THDM-I, and it is almost constant for 
large tan /3. In the model THDM-II, the decay width shows a minimum around tan /3 ~ 7. 
As an example, using the parametrization PI (PII) of THDM-II, the mass m#- = 300 GeV 
and tan /3 = 7, the decay width is obtained T = 5.0 x 10° (3.5 x 10 _1 ) GeV as in Fig. [3] ((4)), 
respectively. 

For the parameter set PI and PII of THDM-I, the branching ratios of charged Higgs 
boson into different decay channels are given in Fig. [HE For a mass value of m A o > 100 
GeV the branching into the channel —> tb becomes dominant as shown in Fig. [6] The 

branching ratios are given in Fig. 1711101 for the parameter set PI and PII of THDM-II. 
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Figure 2. Decay width of H± boson depending on its mass for parameter set PII of model THDM-I. 



Figure 3. Decay width of boson depending on tan /3 for different mass of charged Higgs boson 
for parameter set PI of model THDM-II. 

IV. PRODUCTION CROSS SECTION AT FCC-HH COLLIDER 

The ongoing searches at the LHC rely on specific production and decay mechanism that 
occupy only a part of the complete model parameter space. The cross sections for the 
single production of charged Higgs boson through the process pp —> tH~ + X within the 
THDM-I and THDM-II are presented in Fig. |TT] and [T2j The characteristics of the cross 
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Figure 4. Decay width of H± boson depending on tan /3 for different mass values of charged Higgs 
boson for parameter set PII of model THDM-II. 



M H (GeV) 


Figure 5. Branching ratios to different decay modes of charged Higgs boson depending on tan/3 
and parameter set PI of model THDM-I. 

sections depending on the tan/3 can be seen from these figures as expected from the 
interaction vertices. The cross sections have a minimum around tan/3 ~ 7 for THDM-II 
while it decreases with the increasing values of tan /3 for THDM-I. There exist large values 
of cross section a ~ 37.5 pb for tan/3 = 1, however it is 0.8 pb (1.5 pb) for tan/3 = 7 
within the model THDM-I (THDM-II), respectively. For large values of tan/3 the cross 
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Figure 6 . Branching ratios to different decay modes of charged Higgs boson depending on tan (3 
and parameter set PII of model THDM-I. 
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Figure 7. Branching ratios to different decay modes of charged Higgs boson depending on its mass 
for tan /3 = 1 and parameter set PI of model THDM-II. 

section decreases for THDM-I, while it is in the same level (comparing the cross sections at 
tan (3 = 1 and tan (3 ~ 50) for THDM-II. 

In order to examine the kinematical distributions of associated b(b) quark, the transverse 
momentum distributions of b(b) quarks for the process pp —> ttb(b) + X are presented in Fig. 
fl3l and fl4l for the parameter set PI and PII of THDM-I, respectively. Fig. fT5l and fl6l show 


9 






















M H (GeV) 


Figure 8. Branching ratios to different decay modes of charged Higgs boson depending on its mass 
for tan 0 = 7 and parameter set PI of model THDM-II. 
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Figure 9. Branching ratios to different decay modes of charged Higgs boson depending on its mass 
for tan 0 = 1 for parameter set PII and model THDM-II. 

the pseudorapidity distributions of b(b) quarks for parameter set PI and PII of THDM-I, 
respectively. The invariant mass distributions rritb (t and b(b) quark in the final state) are 
presented in Fig. [T71 and IT81 for parameter PI and PII of THDM-I, respectively. 

Transverse momentum distributions of b(b) quarks are presented in Fig. [T9] and [20] for 
parameter set PI and PII of THDM-II, respectively. Fig. [21] and [22] show the pseudorapidity 
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Figure 10. Branching ratios to different decay modes of charged Higgs boson depending on its mass 
for tan (5 = 7 for parameter set PII and model THDM-II. 
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Figure 11. Cross section for charged Higgs boson production within the THDM-I at FCC-hh 
collider. 

distributions of b(b) quarks for parameter set PI and PII of THDM-II, respectively. The 
invariant mass distributions m t b (top and bottom quark in the final state) are presented in 
Fig. [231 and 1241 for parameter PI and PII of THDM-II, respectively. 

The cross sections for the background processes pp —> ttb(b ) + A", pp —>• ttc(c ) + X and 
pp —> ttj + X are presented in Table U for the center of mass energy y/s = 100 TeV of 
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Figure 12. Cross section for charged Higgs boson production within the THDM-II at FCC-hh 
collider. 



Figure 13. Transverse momentum distributions of b(b) quarks for set PI of THDM-I. 


FCC-I 1 I 1 collider. In the study we consider the background due to one leading 6-jet. Other 
6-jets originating from top quark decays will have different kinematical distributions from 
the top-antitop associated leading one. In addition to the basic kinematical cuts pr > 20 
GeV and \rj\ < 2.5 an invariant mass cut \m t b — m#±| < 0.1m#± is applied for the analysis. 
For the calculation of background cross section A Ub within the invariant mass interval 
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Figure 14. Transverse momentum distributions of 6(6) quarks for parameter PII of THDM-I. 



■nb(b) 


Figure 15. Pseudorapidity distributions of 6(6) quarks for set PI of THDM-I. 

| m t b — m H ±\ < 0.1 m H ± , we assume the efficiency of 6 -tagging to be = 50% and the 
rejection ratios to be 10 % for c(c) quark jets and 1 % for light quark jets since they are 
assumed to be mistagged as 6 -jets. 

In Tables Umvl we present the signal cross sections A as within the invariant mass interval 
| m t b—‘iTiH± | < 0.1 rriH± , for single production of charged Higgs boson in the model framework 
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Figure 16. Pseudorapidity distributions of b(b ) quarks for parameter PII of THDM-I. 



m tb (GeV) 

Figure 17. Invariant mass distributions of b(b) quarks for set PI of THDM-I. 

of THDM-I and THDM-II (for the parametrizations PI and PII) at FCC-hh collider with 
\/s = 100 TeV. 

The final states result from the decays of W bosons for W + W~ + 3 bj et (where we assume 
at least one W boson decays leptonically or both W bosons decay leptonically). In Table 
ED the statistical significances S/\fB (where S is the number of signal events and B is 
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Figure 18. Invariant mass distributions of b(b) quarks for parameter PII of THDM-I. 



Figure 19. Transverse momentum distributions of b{b ) quarks for set PI of THDM-II. 

the number of background events) for the integrated luminosity of L int = 500 fb _1 , and 
different types (THDM-I and THDM-II) and parametrizations (PI and PII) of the model 
are presented. 

In Figures [25] - [26] we present the luminosity requirement for the signal observability 
depending on the mass of charged Higgs boson in the channel - one W boson decays hadron- 
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Figure 20. Transverse momentum distributions of 6(6) quarks for parameter PII of THDM-II. 
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Figure 21. Pseudorapidity distributions of 6(6) quarks for set PI of THDM-II. 


ically while the other decays leptonically (channel - both W bosons decay leptonically) 
for single production of charged Higgs boson within the model framework of THDM-II 
(parametrization PII) at FCC-hh with y/s = 100 TeV. The results are shown in different 
types of lines for different tan /3 values. 
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Figure 22. Pseudorapidity distributions of b(b) quarks for parameter PII of THDM-II. 
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Figure 23. Invariant mass distributions of b(b ) quarks for set PI of THDM-II. 


V. CONCLUSIONS 


Possible extensions of the Higgs sector can be searched for a wide range of parameter 
space in the high energy proton-proton collisions. The ongoing searches at the LHC rely on 
specific production and decay mechanism that occupy only a part of the complete model 
parameter space. The decay modes of the Higgs bosons can be well similar to the background 
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Figure 24. Invariant mass distributions of 6(6) quarks for parameter PII of THDM-II. 

Table I. The cross sections (in pb) for the background processes pp —>• ttb(b) + X, pp —»• ttc(c) + X 
and pp —> ttj + X calculated in the invariant mass range |nitb — rn H ± \ < 0.1 m H ± at the center of 
mass energy yfs = 100 TeV. 


m t b (GeV) 

pp —»• ttb(b) + X 

pp -» ttc(c) + X 

pp —>• ttj + X 

Aa B (pb) 

300 ± 30 

2.99 x 10° 

3.65 x 10° 

2.51 x 10 2 

4.37 x 10° 

500 ± 50 

2.88 x 10° 

3.62 x 10° 

1.23 x 10 2 

3.03 x 10° 

700 ± 70 

1.89 x 10° 

2.40 x 10° 

8.04 x 10 1 

1.99 x 10° 

900 ± 90 

1.34 x 10° 

1.86 x 10° 

5.47 x 10 1 

1.40 x 10° 


Table II. The cross sections for the signal process pp —>• ttb(b ) + X within the THDM-I and 
parametrization PI calculated in the invariant mass range |mtb — | < 0.1m^± at the center of 

mass energy y/s = 100 TeV. 


THDM-I and PI 

A <Js{pb) 

m tb (GeV) 

tan (3 = 1 

tan/3 = 7 

tan /3 = 30 

tan /3 = 50 

300 ± 30 

1.24 x 10 1 

1.30 x 10 -2 

4.09 x 10" 5 

5.02 x 10 -6 

500 ± 50 

2.29 x 10° 

1.97 x 10 -3 

5.02 x 10 -6 

5.96 x 10 -7 

700 ± 70 

5.52 x nr 1 

4.60 x 10" 4 

8.99 x 10" 7 

1.19 x 10- 7 

900 ± 90 

1.13 x 10' 1 

6.59 x 10“ 5 

2.43 x 10" 7 

2.54 x 10“ 8 
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Table III. The same as UH but for the THDM-I and PII. 


THDM-I and PII 

A as{pb) 

m tb (GeV) 

tan p = 1 

tan p = 7 

tan p = 30 

tan p = 50 

300 ± 30 

1.88 x 10 1 

7.16 x 10~ 2 

4.21 x 10“ 4 

7.57 x 10" 5 

500 ± 50 

8.33 x 10° 

5.56 x 10 -2 

1.94 x 10 -4 

5.03 x 10 -5 

700 ± 70 

3.41 x 10° 

1.84 x 10 -2 

1.64 x 10 -4 

3.02 x 10 -5 

900 ± 90 

1.59 x 10° 

1.72 x 10" 2 

9.28 x 10“ 5 

6.32 x 10“ 6 


Table IV. The same as m but for the THDM-II and PI. 


THDM-II and PI 

A (Tsipb) 

m tb (GeV) 

tan p = 1 

tan P = 7 

tan p = 30 

tan p = 50 

300 ± 30 

1.22 x 10 1 

3.57 x 10 -2 

5.22 x 10 _1 

7.24 x 10° 

500 ± 50 

2.45 x 10° 

5.26 x 10" 3 

3.48 x nr 1 

2.02 x 10° 

700 ± 70 

3.90 x 10" 1 

9.37 x 10 -4 

6.08 x 10 -2 

2.97 x 10" 1 

900 ± 90 

1.00 x 10" 1 

1.86 x 10" 4 

1.24 x 10“ 2 

8.59 x 10 -2 


reactions from top and bottom quarks and other sources. If the single production of charged 
Higgs boson associated with top quark is observed at the LHC, one of the following questions 
is to identify the H~tb interaction. The studies on the observables related to the angular 
distribution of charged lepton in the final state and the forward-backward asymmetry can 


be found in 


14] and references therein. Even it seems challenging to measure precisely 


due to the large hadronic background and systematic uncertainties, we look forward to its 


Table V. The same as [U but for the THDM-II and PII. 


THDM-II and PII 

Aa s (pb) 

m tb (GeV) 

tan p = 1 

tan P = 7 

tan p = 30 

tan P = 50 

300 ± 30 

1.70 x 10 1 

2.12 x 10 _1 

3.70 x 10° 

1.12 x 10 1 

500 ± 50 

8.68 x 10° 

6.14 x 10" 2 

1.50 x 10° 

4.96 x 10° 

700 ± 70 

2.59 x 10° 

2.72 x 10" 2 

7.25 x 10” 1 

1.97 x 10° 

900 ± 90 

1.21 x 10° 

1.24 x 10 -2 

5.13 x 10 _1 

1.18 x 10° 
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Table VI. The statistical significances S/^/B at integrated luminosity Lj n j = 500 fb _1 for different 
types (THDM-I and THDM-II) and parametrizations (PI and PII) of the model. The numbers in 
the parenthesis show the results for the channel both TV-boson decay leptonically. 


Model (tan/3 = 1) 

THDM-I 

THDM-II 

m tb (GeV) 

PI 

PII 

PI 

PII 

300 ± 30 

402.65 (230.73) 

612.95 (351.24) 

395.02 (226.35) 

553.09 (316.93) 

500 ± 50 

89.53 (51.30) 

324.52 (185.96) 

95.45 (54.69) 

338.11 (193.75) 

700 ± 70 

26.56 (15.22) 

164.31 (94.16) 

18.76 (10.75) 

124.56 (71.38) 

900 ± 90 

6.50 (3.73) 

91.53 (52.45) 

5.73 (3.61) 

69.53 (39.84) 



Figure 25. The luminosity need to obtain a 3a significance depending on the mass of charged 
Higgs boson and different tan/3 values (solid line for tan/3 = 1, dashed line for tan/3 = 7, dotted 
line for tan/3 = 30, dot-dashed line for tan/3 = 50) within the THDMs. From the decay modes 
ttb(b ) —>• 3bj e t + 2j + l + MET, the final state is accounted for at least 3 6-jets, 2 light jets, single 
charged lepton, and missing transverse momentum. 


exploitation in precision LHC physics and FCC physics scenarios. It is shown that with an 
integrated luminosity of 500 fb _1 at the center of mass energy yfs = 100 TeV of FCC-hh 
collider, the signal can be distinguished from the background for the charged Higgs boson 
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Figure 26. The same as Fig. [25l but for the decay modes ttb(b) —»• 3bj e t + 2l + MET, the final state 
is accounted for at least 3 6-jets, 2 opposite charged leptons and missing transverse momentum. 

mass up to 1 TeV for a large parameter space of two Higgs doublet model. 
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